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We present a numerical study for different discretisations of the static action, concerning cut-off effects and the 
growth of statistical errors with Euclidean time. An error reduction by an order of magnitude can be obtained with 
respect to the Eichten-Hill action, for time separations beyond 1.3 fm, keeping discretization errors small. The 
best actions lead to a big improvement on the precision of the quark mass Mb and Fb^ in the static approximation. 



1. INTRODUCTION 

Euclidean correlation functions in the static ap- 
proximation are known to be very noisy, the noise 
to signal ratio J^ns growing exponentially with 
the time separation. For the Eichten-Hill (EH) 
action P, the law 



together with gauge, cubic and parity invariance 
and locality. Writing the static lattice action 



(2) 



a 



-Bstat -"T-TT, (1) 



noise , . , . 

-Kns = -. 7 oc exp (xo A) , A 

signal 

is roughly fulfilled j2j, with the ground state en- 
ergy of a heavy meson i?stat being linearly diver- 
gent while approaching the continuum limit. 

Here we explore the possibility of reducing the 
exponent in eq. by changing the discretisation 
of the static lattice action. On the other hand we 
want to retain some properties of the Eichten-Hill 
action in order to preserve the same level of 0(a) 
improvement, in particular 

• heavy quark spin symmetry, 

• local conservation of heavy quark flavor 
number, 



wc studied the cases 



Wsix,0) = V 



f(l«'.K 
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Whyp{x,0) = Vhyp 



(3) 

(4) 
(5) 
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where V is the average of the 6 staples around 
the link U{x,0) and Fhyp is the HYP-link 0.2 

2. NUMERICAL RESULTS 

We mainly looked at the gain concerning noise 
reduction and at some scaling properties for 
the new actions. The setup is defined by the 
Schrodinger functional scheme implemented with 
non-perturbatively 0(a) improved Wilson actions 



^For the coefficients ai,a2 and 03 defining the HYP- 
smearing we used the values proposed in 



2 



for the gauge and the hght quark sectors. A more 
detailed discussion of the framework and defini- 
tions of the correlation functions can be found 
in 0]. Our aim is the computation of the Bg- 
meson decay constant as well as other b-physics 
matrix elements. 

2.1. Noise reduction and cut-off effects 

In Fig. n we show i?NS for the static-light axial 
correlator j'^'^.t obtained from an ensemble of 
2500 quenched configurations. This is one of the 
ensembles used for the computation of Fb^ ^toBj, 
in the static limit. The Figure shows that more 
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Figure 1. i?NS for /i''^* from a 24^ x 36, /3 = 6.2 
lattice. See text for the symbols. 

than an order of magnitude can be gained in 
R^s{xo — 1-5 fm) with respect to the EH ac- 
tion (filled circles) by using S'ftat (empty squares), 
which behaves similarly to S'^^t- The picture is 
even better for S^^^ (empty triangles). 

We then studied the scaling behaviour for a set 
of observables. In ^ we reported about the step 
scaling function of the static-light axial current 
renormalisation constant Zj^'^^ (see [5]). Here we 
want to present a study of the quantity 



/rHr/2)l,^o 
/f-*(r/2)|,^^/2 



at ma 



0, 



(6) 



where 9 is the angle defining the periodicity of the 
fermions. Our results are shown in Fig. [2 they 
refer to four different lattice resolutions of an 
volume with L — 1.436ro with vq ~ 0.5 fm. 

Throughout all our computations we defined 
the 0(a) improved correlator /^*''* by using the 
tree level value 1/2 for b^^^^ <6 and the 1-loop 
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Figure 2. Scaling plot for h. Symbols as in Fig.Q] 
(empty circles refer to S^^^). 

results for df^^ . In particular we set 



cT' = -0.08237<?2 + 0(5^), Sst.t = ^™ 



cr'= -O.1164(10)go' + O(go'), ^.tat =5Lt,^> 

stat 



stat 1 
S cA 



--0.090(3)3^ + 0(34), ^,tat = 5, 



HYP 



The first two numbers have been worked out in 
perturbation theory, while the last one has been 
numerically estimated by solving with respect to 
c5j^*(S'^^^) the imphcit equation 



/r(r/2)|^ /i'-*(T/2)|^7^ 



1. 



(7) 



This defines an improvement condition for a dis- 
cretisation (S"^^^), once the improvement co- 
efficients for another discretisation (5'^at) 
known. In Fig. |31 we show our numerical 



results for the coefficients 

stat,(l)/c'HYP 



stat,(l) 



(•^ftat) and 



-A 



{S^^J), taking as input cf-\S^,^,). The 
same data sets used for computing h in eq. © 
have been used here {9 — 0, 6*' = 1/2). The lower 
part of the plot can be regarded as a test of the 
method. 

2.2. The Bg-meson decay constant 

The setup has been used for computing 



■/ J: 



(8) 



for which we have calculated the regularisation 
dependent part of the renormalisation constant 
Zrgi exactly as done in 5 for the EH action. 
Here the i?s-meson boundary to boundary cor- 
relation function /i enters. The proportional- 
ity symbol in eq. ||SJ) summarises volume factors 
coming from the normalisation of the correlation 
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Figure 3. Numerical results for . The 

dashed hne in the lower part is the analytic per- 
turbative result and the error band in the upper 



part is our estimate for c^'^* 



for Ss, 



HYP 



functions, cmp. In addition we introduced 
hydrogen-like wavefunctions on the boundaries of 
the Schrodinger functional in order to minimise 
the overlap with the first excited state; we ar- 
rive at a plateau of length 0.8 fm in <I>rgi and 
the desired ground state matrix element can be 
extracted with confidence. 

The quantity Fb^ ^/f^B^ is related to $rgi via 



0(1/Mb), 



(9) 



Cps(Afb/A^: 

the function Cps (Mb/A-p-j) |S] being known in 
perturbation theory up to and including g^{mt,) 
corrections to the leading order [Jj. 

We computed $rgi for three different lattice 
resolutions {(3 — 6, 6.1 and 6.2) on a L'^ x T topol- 
ogy with T/L = 3/2 and L = 1.5 to 1.9 fm. The 
continuum limit extrapolation of ^rgi, eval- 
uated from S^^^ at xq = T/2, is shown in Fig.^ 
The extrapolated value is 



3/2 



$RGi = 1.74(13). 



3. CONCLUSIONS 



(10) 



We have shown that the problem with the 
statistical precision of correlation functions com- 
puted with the Eichten-Hill action can be over- 
come by changing the lattice discretisation. In 
particular the largest gain is obtained by making 
use of the HYP-links in the static action. Cut-off 
effects for the proposed actions turned out to be 
of the same size as for the EH action, and are 
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Figure 4. Continuum limit of <I'rgi»'o ■ 

in general rather small. The improvements pre- 
sented led to a computation of the quantity $rgi 
in the continuum limit with an error of 7%. We 
are presently reducing this error further by adding 
one more lattice resolution (/? — 6.45). The same 
data can be used for improving the results on 
the renormalisation group invariant b-quark mass 
along the lines described in |H]. Finally, a precise 
determination of Fb^ can be obtained by combin- 
ing the result in eq. H10|l with data around the 
charm quark mass |9I10| . 
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